Reaction, Levy Flights, and Quenched Disorder 
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We consider tlie A + A ^> ib reaction, where the transport of the particles is given by Levy flights 
in a quenched random potential. With a common literature model of the disorder, the random 
potential can only increase the rate of reaction. With a model of the disorder that obeys detailed 
balance, however, the rate of reaction initially increases and then decreases as a function of the 
disorder strength. The physical behavior obtained with this second model is in accord with that for 
reactive turbulent flow, indicating that Levy flight statistics can model aspects of turbulent fluid 
transport. 
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I. INTRODUCTION 

Levy flights have been used to model a variety of physi- 
cal processes such as epidemic spreading [n], self-difFusion 
in micelle systems Q, and transport in heterogeneous 
rocks 0. Levy flights are essentially a generalization of 
ordinary Brownian walks. The normalized step size dis- 
tribution for Levy flights in d dimensions is 



P{r)d'^r 



#■ 



-fdrdn 



(1) 



where r is the step size, / is the step index, Sd = 
27r^/^/((i/2 — 1)!, and rg a lower microscopic step cut- 
off. In the case of / = 2, we recover Brownian motion. 
However, for / < 2, the distribution of step sizes exhibits 
a long-range algebraic tail that corresponds to large but 
infrequent steps, so called rare events. Due to these 
rare, large steps, the mean square step size deviation di- 
verges, and the central limit theorem does not hold. The 
rare, large step events prevail and determine the long 
time behavior. The dynamic exponent z that charac- 
terizes the mean displacement as a function of time by 
{\r{t)\y ~ (const)i^/^ depends on the microscopic step 
index / according to the relationship z = f, indicating 
anomalous enhanced diffusion, that is, superdiffusion. 

It is well-known that quenched random disorder leads 
to sub-diffusive behavior in two-dimensional Brownian 
walks. Levy ffights in such random environments have 
attracted increasing attention recently. The interplay be- 
tween the "built-in" superdiffusive behavior of the Levy 
ffights and the effect of the random environment gen- 
erally leading to sub-diffusive behavior has been exam- 
ined |j, |[ ^. Surprisingly, an e expansion shows that 
in the models of random potential disorder examined 
to date, the dynamic exponent z locks onto the Levy 
ffight index / in any dimension, regardless of the range 
or strength of disorder (notwithstanding some claims re- 
garding divergence- free disorder in [pj). 

The behavior of chemical reactions with random po- 
tential [^, ^ and isotropic turbulence disorder |^, |I^ 
has been examined. Reactants diffusing according to 



Levy flight statistics have also been studied in a model 
of branching and annihilating processes [O . In general, 
potential disorder tends to slow down the diffusing reac- 
tants. Since these reactions typically become transport- 
limited at long times, potential disorder tends to slow 
down the reaction as well. It was found, however, that 
a small amount of potential disorder added to the tur- 
bulent fluid mixing leads to an increased rate of reac- 
tion. This phenomenon of 'superfast' reaction occurs 
because the disorder traps reactants in local potential 
wells, which quadratically increases the local reaction 
rate, while the turbulence rapidly replenishes the react- 
ing species to these regions. As the potential disorder is 
increased, eventually the rate of reaction decreases, due 
to a slowing of the transport. It is interesting to study 
the behavior of reactants following Levy flight statistics 
in quenched random disorder. The question is: Can the 
Levy statistics mimic rapid turbulent transport and so 
lead to superfast reaction? Furthermore, do the reac- 
tions become transport limited or reaction limited at long 
times? 

In this paper, we analyze two different models of the 
disorder. A conventional literature model, which does 
not satisfy detailed balance, is discussed in section II. A 
new model that does satisfy detailed balance is discussed 
in section III. We study the effect of Levy flight statistics 
and quenched random disorder on the simple bimolecular 
recombination reaction in two dimensions. We focus on 
the physically meaningful two-dimensional case because 
the effects disappear above the upper critical dimension 
that is near two and because other, exact methods of 
analysis are probably more appropriate in one dimension. 
Detailed results of the held theoretic renormalization are 
presented. We conclude this paper in section IV. 



II. REACTION IN A COMMON LITERATURE 
MODEL OF DISORDER 



Including the normal diffusion term, the Fokker-Planck 
equation for Levy flights in a quenched potential field has 



been modeled by [Q, ||, |[ 
dG{r, t) 



dt 



D2\7^G{r, t) + Df{^y^G{T, t) 
+V-[G(r,t)Vt/(r)] 



(2) 



where (V^)-* " is interpreted as the inverse Fourier trans- 
form of —k^, which is a spatiaUy non-local integral oper- 
ator reflecting the long-range character of the Levy steps 
with microscopic step index /. Fourier transforms are 
defined as /(k) = j d'^xe^^'^f{'x). The last term on the 
right hand side of Eq. (||) is a drift term due to the motion 
of the walker in the force field. We assume a Gaussian 
distribution of the random potential force field, V{y), 
with correlation 



7 



(nkiMk^)) = ^(27rr<5(ki +k2) 



(3) 



Note that the strength of disorder is parameterized by 7, 
whereas the range is parameterized by y. 
The reaction we are considering is 



A + A 



(4) 



A field theory is derived by identifying a master equation 
and using the coherent state mapping |12(| . The random 
potential is incorporated with the replica trick |0| . The 
action within the field theory is 



d'^x 



dta,{^,t)[dt ~ DfiVy/^ - D2V^ 



-5{t)]a,{^,t) 



A 



d^x / di[2aj(x,t)ai(x,f) 




-Ha-(x,i)a,^(x,i)] - no / (i'^xaj(x,0) 

- J / dhdt2 f (2^)''(5(ki + k2 + k3 + k4) 

^ J Jkik2k3k4 

ki-(ki+k2)k3-(k3+k4) 

"" |ki-l-k2|2+^ 

Xaj(ki,ii)ai(k2,ii)aj(k3,t2)aj(k4,i2) (5) 



Summation is implied over replica indices. The notation 
/i^ stands for /d''k/(27r)'^. 

The concentration of the reactant A at time t, averaged 
over the initial conditions, is given by 



c^(x, i) = lim (a,;(x, i)) 



(6) 



where the average is taken with respect to exp(— S*). 

To apply the field theoretic renormalization procedure, 
the action is recast [^ as 



S 



d^x / dta,(x, t)[dt - ZfDfRiV') 
Jo 

-tif-^Z2D2RDfiiV^ + S{t)]a,{^, t) 
1 



2\//2 



f-''ZxXRDfRfd''xf ' di[2a,(x. 



t)al{-K,t) 





FIG. 1: (a) Diagram representing the propagator. The arrow 
points in the direction of increasing time, and double lines 
represent the bar fields, (b) Normal diffusion vertex D2. (c), 
(d) Reaction vertices A. (e) Disorder vertex 7. 



-fa^(x,i)a^(x,t)] - n'^noR / rf''xa,(x,0) 
-^fi'Z^jRDJR f dhdt2 f (27r)'^ 

^ J Jkikakaki 

xn ,1 ,1 ,1 ski-(ki-|-k2)k3-(k3+k4) 

X j(ki + k2 + k3 + k4) ii^^Tk^y 

Xaj(ki,ii)ai(k2,ti)aj(k3,t2)aj(k4,t2) (7) 

where the renormalization constants Zf,Z2,Z\, and Z^ 
have been introduced to absorb the UV divergences of the 
model. The parameters, DfR, D2R, \r, and 77?, are the 
dimensionless expansion parameters of the model. Since 
the Levy flight term, DfR{\/^)^/^, is the most important 
term, it is chosen as the dimensionless free term, i.e. the 
propagator is [dt — D/i?(V^)'^/^]~^. Note, we are not 
allowed to treat the regular diffusion term as free term, 
as this violates the physics of the scaling and leads to a 
diverging renormalized Dj. Now, the critical dimension 
following from standard power counting of jr, is dc = 
2f + y — 2, and we introduce for an e expansion e — dc — d. 
The scale-setting wavenumber parameter is denoted by 
/i, and we assign the dimensions to the rest of the terms 
accordingly. 

The connections between the renormalized and un- 
renormalized parameters are 



ZfDfR = 


= Df 


flf-^Z2D2RDfR -- 


= D2 


^^f-^ZxXRDfR -- 


= A 


fJ^'^noR = 


= "0 


fi'Z^^RD^R -- 


= 7 



(8) 

To one-loop order, self-energy diagrams and vertex di- 
agrams are summarized in Fig. 1 and Fig. 0. We may 
be tempted to use the momentum-shell renormalization 
procedure here. However, due to the difficulty of regu- 
larizing this action, the first self-energy diagram would 
incorrectly contribute to D2R by the momentum-shell 




FIG. 2: One-loop diagrams: (a) self-energy diagram con- 
tributing to Df, (b) self-energy diagram contributing to D2; 
(c), (d), and (e) vertex diagrams contributing to 7; (f), (g), 
(h), and (i) vertex diagrams contributing to A. 



renormalization, rather than to Dfji by the field theo- 
retic renormalization that is consistent with perturbation 
theory. In the evaluation of the diagram in Fig. 0a, it is 
important to treat the external momentum exactly. If a 
series expansion in the external momentum is performed 
on the integrand, rather than on the result of the inte- 
gral, an incorrect contribution to Z2 arises. Interestingly, 
when the diagram in Fig. 0a is evaluated, only a finite 
contribution to Zf results. Complete calculation shows 



Z2 
Zx 



= 1 



= 1- 



1- 



(2/ - 3)7 



47re 



Z^ = 1 



27re 
IR 
2tt€ 



Xf 



An{f - d) 



(9) 



where a double pole expansion, 1/e and l/(/ — d), is 
used in the calculation of Zx. The use of the double 
pole expansion in two dimensions means we consider / 



just slightly smaller than 2 and y small. As usual, the (3 
functions defined by 

d 

PD2R — tJ-Js-D2R 

give the flow equations in two dimensions as 



dD2R 

dl 
d\R 

dl 
d^R 

dl 



~Pd,r = U- 2)D2r + 



(2/ - 3)jrD2r 



= -/3a„ = (./ - 2)A;, + 



27r 



47r 
47r 



-/3^ 



Zn 



(11) 



where we use the relation /i = A/e', and A is a micro- 
scopic cutoff. Since Zf — 1, the reaction and disorder 
terms do not affect the dynamic exponent, and z = f. 
We determine the long-time decay from the flow equation 
via matching to short-time perturbation theory. The flow 
equations are integrated to a short time such that 

t(r)=iexp[-/ z{l)dl]=tQ (12) 

Jo 

At short times, we find the mean displacement of an 
unreactive particle from {\r{t{l*),l*)\)f = 4£)/t(P), 
and the concentration of reactants from CA{t{l*),l*) — 
[(nofl(r)A'')-i + DfXB.il*)Af-^t{l*)]-\ The long-time 
asymptotic values are given by scaling: (|?'(i)|) = 



3/'*(|r(t(r),r)|)/,andcA(t) 



-21' 



caW*)^*). 



We first investigate the behavior of D2r- As we will 
see below, the fixed point for 7^ is max(0,27re). Using 
this result, we see that D2R flows exponentially to zero 
as long as y < 4(2 — /), when / is near 2. Likely, a higher 
loop calculation would extend the region in which D2R 
flows to zero. Thus, at least within the region in which 
our flow equations apply, D2R always flows to in the 
presence of Levy flights. It is, therefore, unnecessary to 
introduce such a normal diffusion term in the model. 

For e < 0, i.e. in region I of Fig. |^, there is only a 
set of trivial stable flxed points, 7^ = and A|j = 0, for 
the system. The matching procedure gives the normal 
concentration decay as 



CA{t) 



1 



(13) 



In the region of e > 0, 7]^ = 27re is the nontrivial 
fixed point. But depending on the value of <5 = 3/ -I- 
2/ — 6, the matching procedure yields different results. 
For 5 < 0, i.e. in region II of Fig. y, there is no non- 
trivial fixed point for Afl. The corresponding asymptotic 
concentration decay is a little faster than that in region 
I: 



CA{t) 



1 

~o 



1 



A-K\5\DfM- 



(14) 




FIG. 3: Different regions predicted by the flow equations 
for tlie Levy flight system with disorder model I in two di- 
mensions. The flow equation are accurate for small y and / 
slightly less than 2. The flow equations do not apply in region 
IV. 
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4-2f 6-3f 8-4f 

y (fixed f) 

FIG. 4: Decay exponent for the A + A ^r % reaction in the 
Levy flight system with disorder model I in two dimensions: 
CA{t) ~ (const)i"". 

For (5 > 0, i.e. in region III of Fig. §, A|j. = -iirS is the 
fixed point. The asymptotic concentration decay is the 
fastest: 



CA{t) 



1 



AT:5DfKf-HQ \U) 



(15) 



The relationship between the concentration decay expo- 
nent, a, and disorder range, y, at a fixed value of / is 
plotted in Fig H. Note that the strength of the disorder 
has no effect whatsoever on the concentration decay. 

The reader will note that the disorder can never slow 
down the reaction in this model. This is quite an unex- 
pected result, as these reactions are expected to become 
transport-limited at long times, and disorder should slow 
down the transport. Model I, Eq. (0), is somewhat un- 
physical in that this cannot happen due to the lack of 
disorder contribution to Z f. 



III. REACTION IN A MODEL OF DISORDER 
THAT OBEYS DETAILED BALANCE 

Although Eq. (||) is often used in the literature, it 
does not guarantee a long-time Boltzmann distribution 
for G(r, t). That is, this form of the disorder does not sat- 
isfy detailed balance. A more natural form of the Fokker- 
Planck equation for Levy flights in random disorder is 

^^^ = V^-i • [DfVG{v, t) + G(r, t)VV{v)] (16) 

where V-'^"^ is the inverse Fourier transform of — i/c-^^^k. 
Equation (llfl) can be interpreted as a modification of the 
continuity equation to take into account the long-range 
transport induced by the Levy flights; it is derived in the 
appendix. With the same form of the correlation function 
for the potential, Eq. (H), we have 



S = 



Jo 



dt[2a,{x, t)aj{x, t) + d^^{x, t)a^{x, t)] 



-no 



/"d'^xa,(x,0) 



-1 dt.dtj (2^)'*<5(ki + k2 -K k3 

^ J Jkikakaki 

k(-^ki-\, ■ (ki + k2)k3 ■ (kg + k4) 

"" Iki+ksP+i/ 

Xaj(ki,ii)aj(k2,ii)aj(k3,t2)aj(k4,i2) 



k4) 



(17) 



Again, to apply the field theoretic renormalization pro- 
cedure, the action is recast as 



S 



d^^ I dta,{x,t)[dt~ZfDfR{V 





2^|//2 



-(5(i)]a,(x,i) 



J-'^ZxXRDfR / d'^x 



dt 



X [2a, (x, t)a^,{-}i, t) + d'f{^, i)a^(x, t)] 

-fJ-'^noR / d'^xd^{x,0) - -ffZ^jROJ-R / dtidt2 

i2TTf6(ki + k2 + ka + k4) 



kik2k3k4 
k(-^kt^ki ■ (ki -I- k2)k3 ■ (k3 + k4) 

"" |ki+k2P+J/ 

Xaj(ki,ii)ai(k2,ii)aj(k3,t2)a-,(k4,i2) 



(18) 



where, dc — 2 + y , e = dc — d, and the rest of the param- 
eters are the same as in action I. 

The connections between the renormalized and un- 
renormalized parameters are 

ZfDjR = Df 



Zx^RDfR 



A 

no 

7 



(19) 



The diagrams are the same as those in Figs. 1 and 
y, except that the diagrams from normal diffusion, Figs. 
lb and Ep, are not present. A one-loop diagrammatic 
calculation gives 



Zf 



Z^ = 1 



47re 
7r 



+ 



Xr 



27re ■ 47r(/-rf) 

7fl 

27re 



(20) 



where a double-pole expansion of 1/e and l/{f — d) is 
used in the calculation of Zx. 

The dynamic exponent is given by 



f-'-i""'^'^'^ 



(21) 



This suggests that the Levy flights are significantly 
slowed by the presence of the new random disorder term. 
That is, when detailed balance is obeyed, the disorder 
affects the dynamical exponent. 

The P functions give us the flow equations from Zi in 
two dimensions as 



dXR 

dl 
drfR 

dl 



-f3x^ = (/ - 2)\r + 
-PlB = nR 



An 



X2 

in 



(22) 



This action is well-behaved, and there is no regularization 
difficulty. And indeed, both momentum-shell renormal- 
ization and field-theoretic renormalization yield identical 
flow equations to one- loop order. From the flow equa- 
tions, for e > 0, 7i^ flows to oo, indicating that the dis- 
order dominates over the Levy flights. For e < 0, jr 
flows to 0, indicating transport of Levy flights dominates 
over the disorder. This same behavior, that the disorder 
must be adjusted to be compatible with the transport, 
was found in the turbulent reactive flow problem [HOJ. 
An interesting system arises when those two effects are 
competitive, and so we require e = 0, i.e. y = 0. Un- 
der this constraint, jr — 'y/D'j does not flow. Very 
likely, this result holds to all orders. Further analysis 
of the flow equation for Xr indicates two regimes. For 
'y/D'j > 4(2 — f)n/3, or strong disorder, A^ has a stable 
nontrivial fixed point, given by 



A^ = Anif -2) + 



37 
^1 



\*K^-f /Df 



Following the matching procedure, we have 



'-<" ~ A^ {% 



r+1 



(23) 



(24) 
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FIG. 5: Decay exponent for the A -\- A ^ % reaction in the 
Levy flight system that obeys detailed balance in two dimen- 
sions: CA{t) ~ (const)t~'^. The figure is shown for / = 1.9. 
The reaction is transport limited on the solid curve and reac- 
tion limited on the dashed curve. Note that the region a > 1 
corresponds to 'superfast' reaction. 



However, for 7/i?? < 4(2 — /)7r/3, or weak disorder, there 
is no nontrivial fixed point for \r, and we have 



CA[t) 



[4^(2 - /) 



37/D2]A/-27^y 



t 



S + l/a^D') 



(25) 



For the special case of 7/I?? — 4(2 — f)n/3, Xr decays 
marginally, and we have a logarithmic correction to the 
decay: 



CAit) 



31n(t/to) 



87ri:>/A/-2(l + /)t 



2-/ 



(26) 



In the present case, unlike with the action of section 
II, the range and strength of disorder affect the decay 
exponent significantly. The relationship between the de- 
cay exponent and 7 is plotted in Fig. ^j. We see that a 
small amount of potential disorder leads to an increased 
rate of reaction in the Levy flights. But as the potential 
disorder increases further, the rate of reaction eventu- 
ally decreases. This figure is very much similar to the 
one that showed up in the study of reactive turbulent 
fiow M. In fact, if the Levy flight parameters are re- 
lated to the turbulence parameters Q as / = 2 — 5* and 
7 = AnD^^g* these two models predict exactly the same 
concentration decays. As suggested in Q, in order for the 
reaction to occur, multiple reactants must be trapped in 
regions of low potential energy. After the trapped reac- 
tants are depleted, new reactants must be replenished by 
rapid transport to continue the reaction. Certain combi- 
nations of fast transport, / < 2, and disorder, 7, lead to 
'superfast' reaction, a > 1, as shown in Fig. g. Interest- 
ingly, this result means that the inhomogeneous system 
can have a faster reaction rate than the homogeneous, 
well-mixed system. 



IV. CONCLUSION 

We have analyzed the A + A ^ $ reaction in two di- 
mensional Levy flight systems using two models of ran- 
dom disorder. For a common model in the literature, the 
dynamic exponent always locks to the microscopic step 
index /, and the reaction decay exponent varies between 
1 and 2//. This surprisingly unphysical result that the 
disorder cannot slow down the reaction is due to the fact 
that this model does not satisfy detailed balance. For 
a model that does satisfy detailed balance, on the other 
hand, the disorder can and does modify the transport 
properties of the system. When the disorder is adjusted 
to be compatible with the Levy flight statistics, the re- 
action decay exponent first rises above unity and then 
drops to zero as the strength of disorder is increased. 
These results are identical to those from reactive turbu- 
lent flow, and this harmony suggests that Levy flights, 
properly interpreted, can be a viable model of turbulent 
fluid transport. 
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Appendix 



In this appendix, we derive Eq. (16) from a master 
equation. Particle transport by Levy flight statistics can 
be modeled by hopping dynamics of particles in space. 
The hopping rate in the absence of an external potential 
is simply proportional to Eq. (M. In the presence of 
an external potential, we modify this rate in a way that 
satisfies detailed balance Q : 

a(ro ^ r„) - p//roOP(r„-ro)e-'^[^('-")-^('-°)l/2 (27) 

where a(ro — > r^) is the rate at which particles hop from 
To to Tn. The master equation for the Green's function 



of this process is 
dG{r) 



I d'^x a(r - X ^ r)G(r - x) 

- / d'^x a(r ^ r + x)G(r) (28) 



We now expand this master equation to first order in /3, 
i.e. we look for a Fokker-Planck equation that is linear 
in the potential. Noting that P(k) = 1 — r^k-^ , we find 



^--^/fc^G(k)-^/^G(k-h)f(h) 

x[kf + hf -\k-h\f] +0(13'^) (29) 

In real space, we would write this equation as 



dt ^^ ' 2 



(vy^^Gv) 



+G(V2)//V-y(V2)//2G 



(30) 



Since we are interested in how the long-wavelength fea- 
tures of the potential affect the dynamics, i.e. we take 
into account only the small wavevector contributions in 
Eq. (g), we make an expansion of Eq. ( |29| ) for small h. We 
find {kf + hf-\k- h\f)/2 ~ fk^-'^k ■ h/2 for 1< / < 2. 
Realizing that the Green's function should relax to the 
Boltzmann distribution at long times, and that Df and 
/3 may be renormalized a finite amount by a variety of 
irrelevant operators, we conclude that the proper long- 
wavelength theory should be 



dGjk) 
dt 



= -Dfk^Gik) 



-(3Df [ G{k~h)V{h)kf-^k-h (31) 

Jh 

which is exactly Eq. (|l^) with units of /3 inserted. 
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